Regular-to-chaotic tunneling rates using a fictitious integrable system 
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We derive a formula predicting dynamical tunneling rates from regular states to the chaotic sea 
in systems with a mixed phase space. Our approach is based on the introduction of a fictitious 
integrable system that resembles the regular dynamics within the island. For the standard map 
and other kicked systems we find agreement with numerical results for all regular states in a regime 
where resonance- assisted tunneling is not relevant. 
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Tunneling of a quantum particle is one of the cen- 
tral manifestations of quantum mechanics. For simple 
ID systems tunneling under a potential barrier is well 
understood and described, e.g. by using semiclassical 
WKB theory or the instanton approach For higher- 
dimensional systems so-called "dynamical tunneling" [2[ 
occurs between regions which are separated by dynami- 
cally generated barriers. Typically, such systems have a 
mixed phase space in which regions of regular motion and 
irregular dynamics coexist. Tunneling in these systems 
is barely understood as it generically cannot be reduced 
to the instanton or WKB approach. It has been studied 
theoretically @, |, |, |, 0, H, i, El, Ell, El, El, E| and 
experimentally, e.g. in cold atom systems [16|, [l7[ and 
semiconductor nanostructures [18]. A precise knowledge 
of tunneling rates is of current interest for e.g. eigen- 
states affected by flooding of regular islands [19|, [20[, 
emission properties of optical microcavities [2l| and spec- 
tral statistics in systems with a mixed phase space [22|. 

There are different approaches for the prediction of 
tunneling rates depending on the ratio of Planck's con- 
stant h to the size A of the regular island. In the semi- 
classical regime, h <C A, small resonance chains inside 
the island dominate the tunneling process ("resonance- 
assisted tunneling") []j], EH- In contrast, we focus on 
the experimentally relevant regime of large h (while still 
h < A), where small resonance chains are expected to 
have no influence on the tunneling rates. This regime 
has been investigated in Ref. [14], however, the predic- 
tion does not seem to be generally applicable (see be- 
low). Other studies in this regime investigate situations 
[H, [51] , where dynamical tunneling can be described by 
ID tunneling under a barrier, however, in our opinion 
they are non-generic. A generally applicable theoretical 
description of dynamical tunneling rates in systems with 
a mixed phase space is still an open question. 

In this paper we present a new approach to dynamical 
tunneling from a regular island to the chaotic sea. The 
central idea is the use of a fictitious integrable system 
resembling the regular island. This leads to a tunneling 
formula involving properties of this integrable system as 
well as its difference to the mixed system under consid- 
eration. It allows for the prediction of tunneling rates 
from any quantized torus within the regular island. We 



find excellent agreement with numerical data, see Fig. [TJ 
for an example system where tunneling is not affected 
by phase-space structures like cantori at the border of 
the island. The applicability to more general systems is 
demonstrated for the standard map, see Fig. [H 

We consider 2D maps with one major regular island 
embedded in the chaotic region (Fig. [TJ insets) which are 
described quantum mechanically by a unitary operator 
U [24[. Classically the regular and chaotic region are 
separated, however, quantum mechanically they are cou- 
pled. This coupling has consequences for the eigenstates 
of U. While they are mainly regular or chaotic, i.e. con- 
centrated on a torus inside the regular region or spread 
out over the chaotic sea, they do have at least a small 
component in the other region. This is most clearly seen 
for hybrid states (Fig. [2]i). For a wave packet started 
on the ra-th quantized torus (m = 0, 1, ... , m max — 1) 
coupled to an infinite chaotic sea the decay e~ lmt is de- 
scribed by a tunneling rate j m . For systems with a finite 
phase space this exponential decay occurs at most up to 
the Heisenberg time th = ^/A c h, where A c h is the mean 
level spacing of the chaotic states. Introducing purely 
regular states |^reg) and orthogonal chaotic states |^ch) 
the tunneling rate from such a purely regular state can 
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FIG. 1: (color online) Dynamical tunneling rates from a reg- 
ular island to the chaotic sea for the kicked system (2?J : Nu- 
merical results (dots) and prediction following from Eq. J3]) 
(lines) vs l//i e ff for quantum numbers m < 5. The insets 
show Husimi representations of the regular states m = and 
m = 5 at l//i e ff = 50. The prediction of Ref. [l4| for m — 
with a fitted prefactor is shown (dotted line). 
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FIG. 2: (color online) (a-c) The classical phase space cor- 
responding to some quantum maps U, £T reg , and U c h- (d-f) 
Husimi representation of eigenstates of such maps. Eigen- 
states of U have a regular and a chaotic component, as illus- 
trated in the strongest form of a hybrid state (d). Eigenstates 
l^reg) (|^ch)) of U reg (U c h) are purely regular (chaotic). 



be expressed by Fermi's golden rule 7 = (2tt/H) \V\ 2 p cri , 
where p cri = 1/ A c h oc N cri is the chaotic density of states 
and V = (^chl^l^reg) fc> r a time-independent H. For 
a map U one replaces the local average over matrix ele- 
ments in Fermi's golden rule by an average over all 7V c h 
chaotic states and expresses 7 with respect to the time 
period of J7, yielding 

7 = £m 2 > (i) 

ch 

where v = (^chl^l^reg)- The eigenstates of U cannot be 
used for determining the small matrix elements v, as they 
are neither purely regular nor purely chaotic. 

In order to construct purely regular and chaotic states 
we introduce fictitious regular and chaotic quantum maps 
U reg and U cri [25]. Here U reg is regular in the sense that it 

can be written as e _ljF/re s/^ eff , where H reg is a ID Hamil- 
tonian, which is integrable by definition and h e s is the 
effective Planck constant. H reg has to be chosen such 
that its dynamics over one time unit resembles the clas- 
sical motion corresponding to U within the regular island 
as closely as possible (Fig. |2b). The eigenstates l^reg) of 
U reg are localized in the regular region and continue to 
decay into the chaotic sea (Fig. [2^). This is the decisive 
property of l^reg), which is in contrast to those eigen- 
states of U that are predominantly regular but all have 
a small chaotic admixture. The eigenstates |^ c h) of £T ch 
live in the chaotic region of U and decay into the regular 
island (Fig. [2J). 

As |^reg) and |^ch) are eigenstates of different opera- 
tors U reg and C/chj they are not necessarily orthogonal, 
(^chl^reg) = X with < |x| <C 1. In order to apply 
Fermi's golden rule we introduce orthonormalized states 

l^reg) = |^reg),J^ch) = (|^ch) ~ X* l^reg))/ V 1 ~ lx| 2 , 

leading to (^chlVVeg) = 0. We find up to first order in x 
for the coupling matrix element 

V ~ (^chl^ - ^regl^reg), (2) 



which can be inserted into Eq. (pQ). The appearing term 
^2 c h IV ; ch)(V ; ch| is semiclassically equal to the projection 
operator onto the chaotic region. It can be approximated 
as 1 — Preg, where P reg is a projector onto the regular 
island. This yields 

7m « ||(l-P reg )(C/-C/ reg )|C g )|| 2 (3) 

as our main result, which involves properties of the ficti- 
tious regular system U reg and the difference U — U reg . It 
allows for determining tunneling rates from the regular 
state on the ra-th quantized torus to the chaotic sea. 

The most difficult step in applying Eq. (j3]) is the deter- 
mination of the fictitious integrable system U reg: defined 
by a time-independent ID Hamiltonian H reg (p,q). On 
the one hand its dynamics over one time unit should re- 
semble the classical motion corresponding to U within 
the regular island as closely as possible. As a result the 
contour lines of H reg (p,q) in phase space (Fig. [2)3) ap- 
proximate the KAM-curves of the classical map (Fig. [2^). 
On the other hand the function H reg (p, q) should extrap- 
olate sufficiently smoothly to the remaining phase space 
region. This is essential for the quantum eigenstates of 
H reg to have reasonable tunneling tails in the neighbor- 
hood of the regular island. Finding an optimal H reg is 
a difficult task. In fact, it will resemble the dynamics 
within the island with finite accuracy only, due to the 
generic presence of small resonance chains and the com- 
plicated structure of tori at the boundary of a regular 
island. Similar problems appear for the analytic continu- 
ation of a regular torus into complex space due to the ex- 
istence of a so-called natural boundary 0, d, [TO, QjJ . For 
the quantum tunneling problem at not too small /i e ff and 
thus for a finite phase-space resolution, however, such an 
H reg with limited accuracy can be good enough. We will 
discuss below two approaches [26| leading to a sufficiently 
good H reg for the prediction of tunneling rates. Quantiz- 
ing H reg yields the required quantum mechanical opera- 
tor ?7 r eg = e~ lHve z/ heii with corresponding eigenfunctions 
l^reg)- For the numerical evaluation of Eq. (j3]) in Fig. Q] 
it is convenient to replace tTreglVCg) by e - lE m,/Ka \^ g ) 
and approximate P reg ~ Yl l^reg)(^regl> wnere the sum 
extends over m = 0, 1, . . . , \_A/h e & — 1 / 2 J . 

In the following we will discuss the application of 
Eq. (O for ID kicked systems H(p,q,t) = T(p) + 
V(q) J2 n S(t — n), yielding the classical mapping: q t +i — 
q t + T'(p t ), pt+i = Pt ~ V'(q t +i). The corre- 
sponding quantum map over one kick period is U = 
exp[-iV r ($)/fi e ff] exp[-iT(p)/fieff] = U V U T , where ft e ff is 
the ratio of Planck's constant h to the area of a phase- 
space unit cell. We consider a compact phase space with 
periodic boundary conditions for q G [—1/2, 1/2] and 
p G [—1/2, 1/2]. In order to avoid the influence of reso- 
nances and cantori on the tunneling rates we use a system 
containing one regular island with very small resonance 
chains and a narrow transition region to a homogeneous 
chaotic sea. It is obtained by an appropriate choice of the 
functions V'(q) and T'(p) [27[. The phase space is shown 



3 



in the Husimi function insets of Fig. [TJ After determin- 
ing U reg and \ip™ g ) as d esci> ibed in the last paragraph 
we predict tunneling rates by evaluating Eq. (|3]). Fig. [T] 
shows a comparison to tunneling rates, determined nu- 
merically by absorbing boundary conditions at q = ±1/2 
and taking twice the distance between the eigenvalue 
of the m-th regular state and the unit circle. We find 
excellent agreement for the tunneling rates 7 m over 10 
orders of magnitude. The deviations for the smallest 
7 can be attributed to the beginning of the resonance- 
assisted tunneling regime. We determine H reg using the 
Lie-transformation method [28[. With increasing A", the 
tunneling rates following from Eq. (|3]) converge to a con- 
stant value (see Fig. [3fci) and we choose N = 10 for the 
predictions in Fig. [U Note, that for sufficiently high N 
(not shown in Fig. [3^) H reg and the prediction for 7 are 
expected to diverge. 

We now demonstrate that an analytical evaluation of 
Eq. ([3j) is possible for our example system. We define 
functions V(q) and T(p) by a low order Taylor expan- 
sion of V(q) and T(p), respectively, around the cen- 
ter of the regular island [29]. This results in a uni- 
tary operator UyUf with the following properties: (i) 
The corresponding classical dynamics is not necessarily 
regular, (ii) It is close, however, to a regular quantum 
map U reg beyond the border of the island and can there- 
fore be used in Eq. ([3]) instead of U reg . (hi) Within the 
island it has an almost identical classical dynamics as 
U. Therefore (U — UyUf)\ifj r eg) has almost all of its 
weight in the chaotic region and the projection operator 
1 — Preg can be neglected in Eq. ([3j). With the definitions 

1 +£v = e -*s[ v ®-*®] and 1 + e T = e^ 3 ^"^ 
one obtains j m = ||f7y[ey + e T + eye^Uf \i/>™ g ) || 2 - We 
find that typically the third contribution is negligible, 
leading to 



7m «2 / dq\^ eg (q)\ 



+2 



/ dp|V£ g ( 



1 — COS 



1 — COS 



I V(q) ~ V(q) 
' T(p)-T(p) \ 

fteff / 



(4) 



In the last step the sums over the discrete position and 
momentum values have been replaced by integrals, which 
is valid in the semiclassical limit. Agreement with the 
direct evaluation of Eq. (|3]) is found (not shown). If an 
analytical WKB expression for the regular states |^^ g ) 
is known, Eq. (|4]) can be evaluated further. This is the 
case for a different parameter set [30j which yields a tilted 
harmonic oscillator like island embedded in jl chaotic sea. 
We approximate V(q) — V(q) and T(p) — T(p) by linear 
functions and use a WKB ansatz for the regular wave 
function. It turns out that the integral is proportional to 
the square of the modulus of the regular wave function 
at the border of the regular island. We obtain 
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FIG. 3: Predicted tunneling rate, Eq. J3}, normalized by the 
numerical value for m — 0, /i e ff = 1/32 vs order A of H reg 
corresponding to (a) Fig. Q] (A = 10) and (b) Fig. [4] (A = 4). 

as the semiclassical prediction for the tunneling rate of 
the m-th regular state, where a m = (ra + l/2)(A//i e ff) _1 , 
P m — \f\ — a rn , A « 0.28 is the area of the regular is- 
land and c « 1 is h e s independent by a rough semiclas- 
sical estimate. The prediction, Eq. ([5]), gives excellent 
agreement with numerically determined data over 10 or- 
ders of magnitude in 7 (not shown). Let us make the 
following remarks concerning Eq. ([5|): (i) The only in- 
formation about this non-generic island with constant 
rotation number is A/h e s as in Ref. [141 ]. (ii) While the 
term in square brackets semiclassically approaches 1, it is 
relevant for large /i e ff- (hi) In contrast to Eq. (|4|), where 
the chaotic properties are contained in the differences 
V(q) — V(q) and T(p) — T(p), they appear in the prefac- 
tor c via the linear approximation of these differences. 

The paradigmatic model of quantum chaos is the stan- 
dard map {T{p) = p 2 /2, V(q) = -AT/(4tt 2 ) cos(27r<?)), 
which for K = 2.9 has a large generic regular island. 
Absorbing boundary conditions at q = ±1/2 lead to 
strong fluctuations of the numerically determined tun- 
neling rates as a function of h e &, presumably due to can- 
tori. When choosing q = ±1/4, which is closer to the 
island, we find smoothly decaying tunneling rates (dots 
in Fig. [4]). Evaluating Eq. ([3]) gives good agreement with 
these numerical data. Note that this is the first quan- 
titative prediction of regular-to-chaotic tunneling rates 
for the standard map. Here we determine H reg by first 
using the frequency map analysis [31] for characterizing 
the properties of the regular island. This information is 
used to find the optimal 2D Fourier series of order N for 
H reg . The tunneling rates following from Eq. (|3]) show 
for increasing N the expected divergence (see Fig. [SJd). 
For the predictions in Fig. [4] we choose N — 4 as the 
largest order before this divergence. 

We now want to discuss the relation of our approach 
to previous studies. The semiclassical formula presented 
in Ref. [14| (dotted lines in Figs, [lj [4]) deviates from nu- 
merically determined tunneling rates. It works best for 
the case of constant rotation number, which according to 
Ref. [15] is the approximation used in Ref. How- 
ever, it seems to be not generally applicable. The sys- 
tem studied in Ref. [l3[ can be approximated by a ID 
Hamiltonian H reg (p,q) with a cubic potential. Here the 
tunneling path ends far away from the island. In such a 
case our result is also applicable but the use of the WKB 
expression presented in Ref. [13[ is more convenient. In 
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FIG. 4: (color online) Tunneling rates for the standard map 
(K = 2.9) for m < 2. Prediction of Eq. © (lines) and nu- 
merical results (dots) , obtained using an absorbing boundary 
at q = ±1/4 (gray-shaded area of the inset). Prediction of 
Ref. |14|] for m — with a fitted prefactor (dotted line). 

general situations, however, the main contribution comes 
from tunneling to the neighborhood of the regular island 
as seen e.g. from Eq. (jH). We also performed successful 
tests on the tunneling system investigated in Ref. [321 ]. 



In summary, we have derived a quantum mechanical 
formula Eq. (|3]) for the tunneling rates, which involves 
the fictitious integrable system U reg and the difference 
U — U reg . It is the basis for deriving semiclassical expres- 
sions, which we demonstrated with Eqs. (|4]) and ([5]) for 
the case of a fictitious regular system, that is well approx- 
imated by a kicked system. Still there are open questions 
about dynamical tunneling from a regular island to the 
chaotic sea: (i) Which properties of the regular island 
(e.g. size, winding number, shape) and which proper- 
ties of the chaotic sea are relevant in general? (ii) Can 
the approach be combined with the resonance-assisted 
tunneling description and how can cantori be accounted 
for? (iii) How can it be generalized to time-independent 
Hamiltonian systems, in particular billiards? We hope 
that our approach with a fictitious integrable system will 
allow to answer these questions. 
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